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These notes are a summary of our predictions for the new THERA project, related 
to deep inelastic scattering in the region of ultra low x ( x — > 10~ 7 ). We collect here 
^ \ predictions that satisfy two criteria (i) they do not depend on specific features of the 

i 1 model that we have to use to estimate a possible effect; and (ii) they do not contradict 

the HERA data. 
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Abstract 



1 Introduction: Our hopes and main goals at THERA 



1.1 Three domains of QCD at low x 

Deep inelastic scattering is a unique experiment which allows us to take 'snapshots' of the 
constituents inside a hadron at different moments of time with different resolutions. These 
'snapshots' provide the possibility of finding the degrees of freedom (DOF) that are responsi- 
ble for the interaction in QCD and generalize the theoretical approach from the well defined 
domain of perturbative QCD to the unknown non-perturbative ( confinement ) region, where 
the appropriate theoretical methods are still to be determined. DIS allows one to see the 
constituents of size ph 1/Q, where Q is the photon virtuality, at the time t ph 1/mx, where 
x is the Bjorken variable related to the energy (W) of the process (x = Q 2 /W 2 at low x ). 

HERA data as well as theoretical studies suggest that hadrons have qualitatively diverse 
structure in the three different domains ( see Figs, [l] and ^| ): 

1. Perturbative QCD domain where the constituents are of small size and are dis- 
tributed in a hadron with rather low density ( packing factor of these constituents k 
is small ( k < 1 , see Fig. |3| )); 
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Figure 1: Phase dia- 
gram of DIS. k is the 
packing factor which 
is the number of con- 
stituents multiplied by 
its typical area and di- 
vided by the area of a 
hadron. 



2. High parton density QCD domain in which the constituents are still small and we 
can use weak coupling methods, but their density is so large that their packing factor 
k > 1, and so we cannot treat this system of partons using the established pQCD 
methods; 

3. Non perturbative QCD domain in which the QCD coupling is large, the confine- 
ment of quarks and gluons occurs, and new theoretical methods must be developed to 
explore this region. 

Fig. [| illustrates these 'snapshots' of the constituents at different moments of time ( 
different values of x ). One can see three domains with different distributions of the con- 
stituents in the transverse plane. It should be stressed that the distributions do not depend 
on the reference frame, unlike time which differs in different reference frames. 

Each of these domains has its own theoretical problems that can be clarified by THERA 
experiments. The key problems are shown in Fig. [I]. 

1.2 Brief summary of HERA data 

Brief resume of HERA data: these data can be described by models including parton satura- 
tion, but they can also be described without assuming saturation. However, it turns out that 
all predictions of asymptotic hdQCD have already been seen in HERA data. This fact is so 
impressive and convincing that we, personally, think that HERA has reached a new regime 
of high density QCD However, the situation is still non-conclusive as is illustrated by 
Fig. |] which shows the value of the saturation scale in HERA and THERA kinematic region. 
One can see that Q s (x) < 1 GeV for HERA. This low Q s (x) indicates that HERA data can 
be described by other approaches without saturation, for example, by some models that 
include a smooth matching between "soft" and "hard" interactions. However, at THERA 
Q s {x) is larger and the hdQCD interpretation of the data will be cleaner. We illustrate this 
point with two figures that follow. 
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Figure 2: 'Snapshots' of the 
constituents of different size 
(r « 1/Q) at different val- 
ues of x. Large circles in- 
dicate the hadron and small 
ones the constituents. 
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1.3 Main idea 

As we have discussed, we face two challenging problems in the region of low x and low Q 2 
which is now being investigated at HERA : 

1. The matching of "hard" processes, which can be successfully described using per- 
turbative QCD (pQCD), and "soft" processes, which should be described using non- 
perturbative QCD (npQCD); 

2. Theoretical description of high density QCD (hdQCD). In this kinematic region we 
expect that the typical distances will be small, but the parton density will be so large 
that a new non perturbative approach needs to be developed for dealing with this 
system. 

The main physical idea, on which our approach is based is |J: 

The above two problems are correlated and the system of partons always 
passes through the stage of hdQCD ( at shorter distances ) before it proceeds 
to non-perturbative QCD and which, in practice, we describe using Reggeon 
phenomenology. 

1.4 Status of theory 

Parton saturation as well as other collective phenomena typical of the high parton density 
system, is not an additional postulate of QCD, but follows from the QCD evolution equations 
in the kinematic region associated with high parton density. Therefore, it is very important 
to have a clear understanding what can be proven theoretically. 

In DIS at low x, one can find a system of high density partons, which is a non-perturbative 
system due to high density of partons, although the running QCD coupling constant is still 
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Figure 3: The parton pack- 
ing factor k as function 
of Q 2 and Bjorken x in 
the GRV'98 parameteriza- 
tion of the solution to the 
DGLAP evolution equation. 
The GRV'98 parameteriza- 
tion describes all available 
data from HERA . 
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small ( as(r±) C 1 ), Such a unique system can be treated theoretically 
stressed that the theory of hdQCD is now in very good shape. 
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Figure 4: The estimates 
for the saturation scale for 
HERA and THERA kine- 
matic region. 
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Two approaches have been developed for hdQCD. The first one is based on pQCD 
( see GLR and Mueller and Qiu papers in Ref. ) and on dipole degrees of freedom ||]. 
This approach gives a natural description of the parton cascade in the kinematic region for 
k < 1 and up to the transition region with k ~ 1 ( see Fig. D). 

The second method || uses the effective Lagrangian suggested by McLerran and Venu- 
gopalan , this is a natural framework to describe data in the deep saturation region, where 
k 3> 1 (see Fig. 0). As a result of intensive work using these two approaches the non-linear 
evolution equation which has the following form has been derived 
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where a el (rj_, b t , x) is the elastic scattering amplitude for a dipole of size r± at energy ex 1/x 
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The pictorial form 



of Eq. ([I]) is given in Fig. |5] which shows that the physics underlying this equation has a 
simple meaning: the dipole of size xiq decays in two dipoles of sizes xu and 202 • These two 
dipoles interact with the target. The non- linear term which takes into account the Glauber 
corrections for such an interaction, Eq. (|l|) is the same as the GLR -equation [0] but here 
it is given in the coordinate representation. The coefficient in front of the non-linear term 
coincides in the double log with the one calculated in Ref. 0. 
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Figure 5: The picto- 
rial form of the non- 
linear evolution equa- 
tion in the hdQCD 
kinematic region. 



We wish to stress that this equation which includes the Glauber rescatterings , has definite 
initial conditions and has been derived by both methods (see Refs. @, The model, which 
we will describe in the next section, provides both the correct initial conditions for Eq. ([!]) 
and also serves as a good first iteration. This iteration reproduces the main features of the 
solution, and it is only necessary to repeat the iteration procedure two or three times to 
obtain a correct solution for x < 10~ 7 . 



2 Our model 



2.1 General description 

As was shown in Ref. |J, the correct initial condition for Eq. (P is actually the Glauber- 
Mueller formula [10[ [U], [12|] for the rescattering of the colour dipole, namely 



o-di P oie(r±,x) 



d 2 b t Ima el (r±,x; b t ) 
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where 

a el (r ± ,x;b t ) =i 1 - e — , (3) 



The opacity 0(rj_, x; 6 t ) is defined as 



fi(r ± ,x;6 t ) = -^xG DGi ^(x,-;6 t ) , (4) 
3ttR z 



where G DGLAP (x, 4; M = G DGLAP (x, 4) • SY&t) and G DGLAP (x, A) is the solution of the 
linear DGLAP evolution equation, and S(bt) is the profile function for the impact parameter 
distribution of the gluons in the target. The origin of this function is non-perturbative, and 
it is normalized in Eq. ([|) by the condition S(b t = 0) = 1. 

For the solution of Eq. (|]) one should fix the value of initial x = Xq (y = yo) and use 
a el (r±,x = x ;b t ) as a starting iteration of Eq. (lip. In our model we suggest a different 
approach, namely we use Eq. (0) and Eq. (||) as the first iteration of the Eq. ([[]) including 
their x-dependence. Therefore, the result of the second iteration of Eq. ([!]) can be written 
in the form: 



4(r ± ,y = Hl/x);b t ) = ^ r\ f dy' f ^ { 2 af (r' ± , y'; b t ) - (af(r' ± ,y';b t )) 2 } 

7T J Jr T 



(5) 

where we assumed that r' ± 3> r±. This assumption corresponds to the Leading Log Ap- 
proximation of pertutbative QCD LLA, which has been used in the derivation of Eq. (Q) and 
Eq. (g). Substituting Eq. (§ in Eq. (g) we obtain 

al l (r ± ,y = Hl/x);b t ) = r\ J dy' ^ ( 1 - e i J , (6) 

where Q G (r ± , y- b t ) = 2Q(r ± , y; b t ) of Eq. (§ Q. 

The physical meaning of Eq. (|6]) is transparent. The dipole of size r± decays into two 
dipoles which interact with the target. Eq. (|6|) describes the rescatterings of these two 
dipoles. On the other hand, in pQCD this state is the qqG state. Since we assume the size 
of qq system to be much smaller than the size of the dipoles in the qqG state, Eq. (|) relates 
to the passage of the gluon through the target. Using the relation between the gluon-target 
cross section^] and the gluon distribution 

A 2 

a G = 2 J d\a e 2 l (r ± ,y = \n(l/x);b t ) = ^a s (Q 2 )xG(x,Q 2 ) 



one obtains the Glauber-Mueller formula for the gluon distribution [12 

^(«,(« = A [ ld 4 l *<L /^(i-e-f). (7) 
7r 4 Jx x' Ja/q 2 r'l J V / 

Eq. (0) and Eq. (0) are the main formulae that we use in our estimates of the collective 
phenomena in DIS. 

Actually, fl G /n = 2iV c 2 /(iV c 2 - 1) = 9/4(7V c = 3) _> 2(N C > 1). 

2 In principle, the gluon - target cross sections can be measured using the graviton as a colouless probe. 
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2.2 Advantages and disadvantages of the model. 



The main advantages of our model follow directly from the way it has been constructed. Our 
model reproduces the DGLAP limit for r\ < r 2 saturation « 1/Q 2 S , gives a good approximation 
to the solution of Eq. P) for x > 10 -6 , and it preserves the relation between elastic, quasi- 
elastic ( diffraction) scattering and multi particle production in DIS based on the AGK 
cutting rules [|13|] (see Ref. U for details). 

The main problem relating to our model is the fact that the evolution equation Eq. ([]]) 
has only been proven in the leading ln(l/x) approximation of pQCD where we consider 
as ln(l/x) ~ 1 while as <C 1- This approximation does not insure the accuracy of calculation 
for present accessible energies. On the other had, our model cannot be correct at low x and 
it is only suitable to describe DIS for x > 10~ 6 , where the model gives the second iteration 
of Eq. (|l|). For smaller values of x we require higher iterations. 

In Ref. |T^] we have already shown that our model gives a good approximation of Eq. ([[]) . 
Fig. H illustrates this point. Hence we can safely use our model for estimates of the collective 
phenomena even in THERA kinematic region. 
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Figure 6: Figs.6-a - 6-b show the calulation for N = Ima el at b t = (o-di po i e = 
J d 2 b t N(x,r±;b t )) in our model ( full curve ) and the solution to nonlinear equation (see 
Eq. §)). 

Table 1 provides a guide for the different processes which we described in our model for 
kinematical range at HERA. 

2.3 Phenomenological parameters of the model: 

Before discussing the applications at THERA we list the parameters that we use to fit the 
data at HERA. 
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Reaction 



Table 1 

Q 2 (GeV 2 ) x References 



F 2 (x,Q 2 ) 

xG{Q 2 ,x) 

dF 2 /dlnQ 2 



0^65 
1 ^65 
1 ^ 65 
1 ^ 65 



< o.oi ra 



<0.01 [0,0 

< o.oi rmi 



< 0.01 IS, 23 



diSS 



atotdr) 0;0-20 < 0.01 H 

<r% f 5 ^ 65 < 0.01 |g 



1^65 <0.01 21 



a(Yp ^ J/^+p) 4- 65 < 0.01 [0, 123 



slope B(Yp -> J/tf + p) 4- 65 
slope -B(7*p — > p + p) 5 4 65 



< 0.01 22 



< 0.01 22 



2.3.1 R 2 - size of the target. 

The size of the target enters the impact parameter profile of the target which we take in the 
Gaussian form: 

S(b t ) = -L c -£ . (8) 

The HERA data for photo production of the J/\l/ - meson as well as CDF data on double 
parton cross section, leads to the value of R 2 = 5 4- 10 GeV~ 2 . R 2 is a parameter fiied to 
describe of the experimental data. Note, that the value of R 2 = 8.5Ge^~ 2 was taken for all 
reactions that we have described. 



2.3.2 Qq = l/r 2 ep - separation parameter. 

As we have discussed we can only rely on our model for the saturation effect ( see Eq. (0) 
- Eq. (|D ) at rather small distances (r± < r s ^ p ) or , in other words, at large virtualities 
of the incoming photon Q 2 > Qq. We have commented on the value of r^ ep , but in practice 
we used Qq = 0.6 4 1 GeV 2 and tried to estimate how our fit depends on the value of Qq. 
Therefore, the result of our calculations should be read , as "the shadowing corrections from 
short distances r± < 1/Qq gives this or that 



2.3.3 Solution of the DGLAP evolution equations. 

We attempted to use all available parameterization of the solution of the DGLAP evolution 
equations [|4], p5|j , but we prefer the GRV parameterization [p6[1 . The reason for this is very 
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simple: the theoretical formulae, that are the basis of our model, were derived in double log 
approximation of pQCD, and the GRV parameterization is the closest one to the DLA. 



3 Predictions for THERA 

3.1 The unitarity bound in THERA kinematic region 

We start from the prediction which, in principle, does not depend on the exact form of the 
correct evolution equation and/or on the particular model, namely, from unitarity bound 
for F 2 and xG(x, Q 2 )E7\. This bound stems from a simple formula for the DIS cross section 



m m 

<r(r/*P) = J o dz J d 2 r ± \^(z,r ± ;Q 2 )\ 2 (Tdi P oie(xB,rl) , (9) 

where (Jdipoie(xB, is the total cross section of the qq -dipole of size r± with the target; 
\1/ is the wave function of the qq -dipole in the virtual photon. This wave function is well 
known ||T2| , |16| and for transverse polarised photon \^t{z, r_i_; Q 2 )\ 2 is equal 



„,em A/" 

\y T (z,r ± ;Q 2 )\ 2 = xE^^ + ll-^Q 2 *!), (10) 

where K\ is the modified Bessel function, Q 2 = z(l — z) Q 2 , Nf is the number of massless 
quarks and Zf is the fraction of the charge carried by the quark. 

It was shown in Ref. |27j that in the DGLAP limit the essential r± in Eq. (H) are 
larger than 2/Q ( > 2/Q) and the integral over z can be taken, namely, 

Finally, using the relation between the total cross section and F 2 structure function 

4 7r 2 a em 

0{-fp) = —^z—F 2 (x B ,Q 2 ) (12) 

one obtains 

N f ,oo AJl AT N f 



AT J roo A r ^ AT - ' r 

F 2 (xb, Q 2 ) = ^ £ q /, = E d\lmai pole (x B ,rr,b t ) 

" 1 ± 1 



(13) 

Taking the derivative with respect to InQ 2 and using the weak form of the unitarity 
constraint (Im a d \ vole {x b , t j_; b t ) < 1 ) we obtain 

dF 2 (x,Q 2 ) 

dlnQ 2 ~ 3tt 2 ' 1 ' 

R 2 in Eq. fll"4]) is the region of convergence for the integral over bt in Eq. flU]). In principle, 
R 2 grows with x but model estimates |2^| as well as experimental data |28| show only mild 
x dependence. 
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Fig. |^ shows that the F 2 slope approaches the unitarity bound at least in the GRV'98 pa- 
rameterization of the solution to the DGLAP evolution equation. A violation of the unitarity 
bound in the DGLAP equation indicates that the shadowing corrections are unavoidable in 
the THERA kinematic region. The real size of these corrections is much larger than we can 
see from the violation of the unitarity bound, since the system starts becoming dense at 
densities lower than that which follows from the unitarity constraints. In other words, shad- 
owing corrections lead to a considerable suppression in the deep inelastic structure function, 
at densities lower than originates from the unitarity constraints. 



3.2 The unitarity bound for DIS with nuclei. 

Eq. (jlj]) looks better for DIS with nuclei, since the radius of the nucleus is large and a 
shrinkage of the diffraction peak induced by SC will be very small and we can neglect it. We 
plot the F 2 a = AFffi LAP and unitarity bound in Fig. || 

One can see that for DIS with nuclei we should see the collective phenomena at x ~ 10~ 4 
and at rather large value of Q 2 ~ 3 4- 5 GeV 2 . 



3.3 Scaling violation in the F2 - slope 

The careful analysis of the HERA data on the F 2 - slope ( 9F q\^q2 ' ) > given in Ref.fll8l, shows 
that (i) our saturation model as well as other models of this type ( see [2S] for example 
) , are able to describe all experimental data ; and (ii) such a description cannot be very 
conclusive since other approaches are equally successful. The saturation hypothesis leads to 
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Figure 8: The unitar- 
ity bound 
and the DGLAP pre- 
dictions for the F 2 a - 
slope at different val- 
ues of x and atomic 
number A. 



&F q^q2 ] oc Q 2 R 2 for Q 2 < Q 2 (x). The HERA data ||, |lj show such behaviour, but we 



cannot distingish this saturation behaviour from the vanishing of F2 on the soft scale, which 
follows from the fact that the total photoproduction cross section is finite at Q 2 — > 0. There 
are two reasons for this uncertainty: (i) the soft scale is not so soft and typical transverse 



momentum in the soft Pomeron could be as large as 2 GeV ||32|| ; and (ii) the saturation scale 
is rather small Q 2 {x) = 1 4- 2 GeV 2 in HERA kinematic region. 

One can see from Fig. || that THERA will allow us to distinguish between a mixture of 
soft and hard Pomerons (DL curve in Fig. |S] ) and our model for gluon saturation ( GLMN 
curve in Fig. || ). The difference between the DGLAP approach ( CTEQ5 curve in Fig. 
|S] ) and our predictions is concentrated in the region of small Q 2 ~ 1 4- 2 GeV 2 , but we 
recall that the corrections to the CTEQ5 parametrization due to high parton density effects 
reaches the value of about 30-40% at THERA energies. These estimates are an alternative 
way of saying, that at THERA energies we expect large SC theoretically, and a DGLAP 
approach can absorb these corrections in the initial nonperturbative gluon distributions at 
HERA energies, but it would be a more difficult task in THERA kinematic region. 

3.4 Energy dependence of J/^ production. 



It was shown in [[3^, that the energy behaviour of the J/^ production is very sensitive 

to the value of the shadowing corrections. It turns out that the large uncertainties due to our 
poor knowledge of the wave function of vector mesons contribute mostly in the normalization 
of the cross section, while the energy slope is still a source of the information on SC. In Refs. 



34 1 we showed that the shadowing corrections provide a natural explanation of the 
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;ure 9: Predictions for the different parameterizations at fixed low values of x. 
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Figure 10: A compari- 
son between the energy 
dependence of differ- 
ent models (see text 
) for the integrated 
cross section of J/ty 
photoproduction. The 
available experimental 
data points are con- 
fined within the inner 
window. 



experimental energy behaviour for J/\l/ photo and deep inelastic production. However, the 
available data do not enable us to exclude explanations based on the mixture of soft and hard 
Pomerons, and the idea that this process is hard . Fig. [10] gives our prediction for THERA 
kinematic region for three approaches: the SC calculations in our model (GRV98SC), the 
soft + hard Pomeron model (DL2P)|35| and the DGLAP approach based on the CTEQ 
parametrization (CTEQ5NSC). 

3.5 Shrinkage of diffraction peak for production of vector mesons. 



The experimentally observed shrinkage of the diffraction peak in photoproduction of J/\l/[ 
is a direct indication that this process is not a simple hard process that can be described 
in the DGLAP approach . Indeed, one of the most well established properties of the hard 
processes is the fact that the t-dependence is independent of energy (x) (see for example 
Ref. ||22|| . There are two possible explanations: (i) the first one is the SC which lead to 



x-dependence of the £-slope[B2] and (ii) the second is based on the contamination of the 
j/\[r production by the soft processes for which the shrinkage of the diffraction peak is a 
phenomenon that is well established both theoretically and experimentally. 

It should be stressed that the above two approaches have different predictions for the 
x-dependence: SC lead to the t-slope which increases at higher energies ( lower x), since the 
contribution of SC grows with energy. For the mixture of soft and hard processes, the role of 
soft ones diminishes at higher energies, and as a result the value of effective a' e ^ decreases 



In Fig. 11 one can see this effect for our model. This figure also shows that we cannot 
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Figure 11: The energy 
(x) dependence of the 
for- 
ward differential slope 
of J/ty photoproduc- 
tion.ZEUS data ^ 
and our model calcu- 
lation with several val- 
ues of R 2 . From 
this picture we chose 
the value of R 2 = 
8.5 GeV' 2 for the typ- 
ical proton size in our 
model. 



describe the ZEUS data regarding the value of the i-slope. The reason for this may be due 
to our under estimating the value of SC in our model. However, one lesson we can learn from 
Fig. [II]: THERA will clarify the question which mechanism works. The important thing to 
emphasize once more is that the measurement of the shrinkage of the t-slope will provide 
reliable information on the deviation from the simple DGLAP approach.lt is especially im- 
portant to observe such shrinkage in the DIS diffraction production of J/\l/ and other vector 
mesons. 



3.6 Maxima in ratios 

Preparing this paper we tried to find improved observables which will be sensitive to the 
saturation scale. We study the Q 2 behaviour of the ratios:^/^ and F^fF® for longitudinal 



and transverse structure function for inclusive DIS and for diffraction in DIS Pq| . In Fig. [12 
some examples of these ratios are plotted. We found that these ratios have maxima at 
Q 2 = Qmaxix) which increases with x as x — > 0. It appears that Q max (x) is a simple 
function of the saturation scale Q s (x). In the THERA kinematic region Q 2 max [x) is large 
Q 2 max {x) ~ 6 -T- 7 GeV 2 . Such a large value of Q 2 max {x) makes our calculation more reliable, 
and we expect that the measurement of this maxima in the THERA kinematic region will 
enable us to extract the value of the saturation scale from the experimental data. 
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Figure 12: The ratio of F® / F<f as function of Q 2 at different values of x and the behavior 
°fQliax( x ) as function of x. 



3.7 Higher twist contribution. 

One of the most challenging problem of QCD is to understand the higher twist contributions. 
The present approach to DIS is based on two main ideas: (i) the DGLAP evolution equation 
for leading twist contributions and (ii) the firm belief that higher twist contributions are 
small in the whole kinematic region, when we start QCD evolution from the large value of 
Q 2 = Ql ~ 1 — A GeV 2 . In recent years it has been proven that there is no ground for such 



an assumption. It was found |39| that the anomalous dimension for the higher twists is much 
larger than for the leading one in the region of low x. It turns out that if we write the deep 
inelastic structure function in the form 



F 2 (x,Q' 



M 2 

%Q 2 ) + ^F 2 HT i 



x,Q 2 



F 2 HT i 



x,Q 2 ) 



F 2 LT (x,Q 2 ) x xG(x, Q 2 ) at x — > 0. Therefore, the experimental observation 



of the higher twist contribution, is one of the most challenging and important problems in 
DIS, as well as in QCD at large. The attractive feature of our model is the fact that it leads 
to higher twist contributions in accord with known theoretical information. Our calculations 
confirm the result of Ref. |3(J that there is almost a full cancellation of the higher twist 
contributions in F 2 in spite of the fact that they give substantial contributions separately to 
Fl and Ft as well as to F D . In the THERA kinematic region ( at x ~ 10 -5 ) we expect the 
higher twist contributions to be of the same order as the leading twist at sufficiently high 
value of Q 2 (see Fig. O) This high value of Q 2 insures us that our calculations are reliable. 
Therefore, we believe that THERA has a good chance to measure higher twist terms and a 
new era of DIS will open, that will include a systematic study of higher twist contributions. 
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Figure 13: Different twist contributions to the various structure functions for DIS on the 
proton: leading twist (at high Q 2 ) - dashed line, next-to-leading - dotted one, exact structure 
function - solid curve. 
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4 Resume 



We presented here our estimates for the possible manifestation of saturation in the THERA 
kinematic region. We believe that HERA has reached a new QCD regime: the high parton 
density QCD domain where incorporating new collective phenomena is essential for 
understanding it's physics. We argue here that data from THERA will be able to show 
that we have reached this new regime, and will allow a systematic study of the QCD parton 
system with large parton density. 

We hope that our estimates will help to plan the experimental strategy for the THERA 
project. 
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